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1. Concept oF Dyapics AND POLYADICS. 


A sysTEM of Vector Algebra in N dimensions following Gibbs is 
based on a set of N unit vectors @), -- -, such that the dot products 
€;-@, are unity when subscripts are equal, otherwise zero. Any other 
vector @ may be expressed in terms of the fundamental units, 


& = a0; + + --- + ayey. (1) 


Two vectors a and b written together with no dot or other sign between 
them constitute a dyad ab called the indeterminate product of a into b. 
A dyad e,e, will be called a fundamental dyad or dyad unit. If we 
have any dyad ab the first vector a is called the antecedent, the second 
the consequent. Both a and b will be called factors of the dyad. A 
dyad is zero when and only when all products of scalar elements 
a,b; are zero, that is, when one of the factors vanishes in all its ele- 


ments. 
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A dyadic is a sum of dyads, or, with no increase in generality, a sum 
of dyads each multiplied by a scalar factor. The elements of a dyadic 
are the sums 2a,b; of corresponding products taken from each dyad 
term, where, for an element A; of a dyadic A the subscripts 7 and k 
are constant and the summation is over the various dyad terms. 
Thus if 

A = a,b; + abe + + (2) 
we have 
= + + + andar. (3) 


It is frequently of use to write the dyad ab as symbolical of the dyadic 
and a,b; as symbolical of the element A;, omitting the subscript h 
and the sign of summation with respect to h.1 

A dyadic may be written in terms of the dyad units e;e; se 


A = (4) 


where both subscripts run from 1 to N. 

In a similar manner several vectors abe: --g written with no sign 
between them constitute a polyad. When we wish to indicate that 
the polyad is of order K, that is, it is the indeterminate product of 
K vectors, we may call it a K-ad. The vectors a, b, ¢, ete., in order, 
will be called the first, second, third, ete. factors of the polyad. A 
polyadic is a sum of polyads, all of the same order. The polyads 
e,e,e,::-e, to K factors will be called the fundamental polyads of 
order K or fundamental K-ads. It is evident that any polyadic may 
be expanded in terms of the fundamental polyads of its own order. 
The elements of a polyadic are the scalar coefficients in this expansion. 
Thus if we have K subscripts p, q,---, 8, all of which run from 1 to N, 
we may write a K-adic A as 

A = . (5) 
It is thus apparent that in general a K-adic depends on N*¥ scalar 
elements. These may be primarily regarded as forming a K-dimen- 
sional block, and most conveniently written by means of adjacent 


squares when special values have to be assigned. For example if 
N = 2 and K = 4 we should have the scheme of elements 


Aun, Aji Ajon, 


Ane Ajon1, 
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as one possible way of representing a four-dimensional assemblage of 
elements on a flat surface. 

A polyadic is defined to be zero when and only when all its elements 
are zero. Since our fundamental conception of a polyadic is as a sum 
of polyads rather than as a mere aggregation of scalars, the equation 
A = Ois to be thought of as equivalent to N* equations of the form 


It is frequently useful to indicate such a set of equations symbolically 
as ab--g = 0, where, as before, the polyad is written symbolically for 
the polyadic. ! 


2. Potyapics as VECTORS IN SPACE OF HIGHER DIMENSIONS. 


Investigations on polyadics may be distinguished according to 
whether the K-dimensional character is important or not. If we agree 
on some definite order among the fundamental polyads we may write 


E.= (8) 


and let r run from 1 to N¥ while each of the K subscripts on the right 
runs from 1 to N. With A,= A,,.., we shall then have 


A= + (9) 


where, for convenience, has been put for N¥. 
The polyadic A thus takes the form of a vector in space of n dimen- 
sions. This concept is justified if we introduce the multiple dot 


product defined by 


where the colon in every case indicates K-tuple dot product. In 
words 

Definition. The K-tuple dot product of two K-ads is the product 
of dot products of corresponding factors. 

The K-tuple dot product of two K-adics is the sum of K-tuple dot 
products term by term. By virtue of the distributive law of multipli- 
cation the relation (10) is symbolical of the K-tuple dot multiplication 
of two K-adics. When K = 2 we have Gibbs’ double dot multiplica- 
tion. 

By (8) and (10) it is evident that E;: E;, is unity when 7 = fk, other- 
wise zero. Thus the fundamental polyads behave with respect to 
multiple dot product as do unit vectors @1, @,-:-@, in space of n 
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dimensions with respect to ordinary dot product. Therefore all 
properties of such vectors which pertain to dot multiplication alone 
will go over into properties of polyadics and multiple dot product. 

To say the same thing in another way, the fundamental polyads 
€,@,°:e, form a normal orthogonal system with respect to multiple 
dot multiplication, in the sense that the product of two unlike polyads 
of order K is zero while the product of a polyad by itself is unity. 
The notation E,. for the polvad is convenient but not essential. The 
analogy with e,, while of much value in a formal sense, is temporary 
in character, to be employed or laid aside according as we wish to 
forget or to emphasize the A-tuple nature of the polyads. 


3. THe FUNDAMENTAL IDENTITY OF Dot MULTIPLICATION. 


Let there now be a set of N¥ or n polyadics of order K as Aj, Ae,-:: 
A,, arbitrarily chosen. Any one of these, as A;, may be expanded in 
the manner of (9), thus 


A; + + + QinEn. (11) 


Let there also be another polyadic which we may on occasion call 
A,..; but which it will be more convenient at present to call M and to 
expand in the form 


M = M:E,+ M.E,+---+ M,E,. (12) 


If we now take the multiple dot product M:A;, remembering the 
principle of orthogonality pointed out in the last article, we shall have 


M:A;= Mrain. (13) 


There are n equations of this form, one for each of the polyadics 
A,---A,. Together with (12) we thus have n + 1 equations linear 
in the scalars M;---M,. It is true that (12) is a K-adic equation 
and is itself equivalent to n scalar equations. But since all the 
equations (13) are merely scalar equations, it is not hard to see that 
the determinant 


M: Ap, U2}, Gen (14) 


M: A,, Ann 


must vanish; because when polyadics have to be multiplied by scalars 
only, all the laws of ordinary algebra are obeyed. 
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We have now to develop this determinant by the elements of the 
first column, yielding the identity 


CoM + C,M: A, + C.M: A, +---+C,M:A, = 0 (15) 


where the scalar Cy and the polyadics C,---C, are the cofactors of the . 


elements into which they are multiplied. Evidently Co is the determi- 
nant of the n? elements a;,. To exhibit the character of C;, Co, etc., 
let c;, denote the cofactor of a;, from the n-rowed determinant Cp. 
On developing (14) it appears that 


+ +: ¢inEn) (16) 


whence the ©; are polyadics whose elements are the negatives of the 
n — 1 rowed cofactors from C5. 

So far the n+ 1 polyadics M, Aj, As,---,A, have been quite arbi- 
trary. If, however, the A; are linearly independent, their determinant 
Cois not zero and we may introduce the important new set of polyadics 
A’, A’, defined by 


CoA’; = —C; (17) 
or by (16) 


0 


The set A’; is said to be reciprocal to the set Aj. 

By multiplying the expansions (11) and (18), remembering that 
cx is the cofactor of a;, from the determinant Co, we have the relation 
A’;: A; = 1 when subscripts are equal, otherwise zero, which we may 
also express by saying that a set of n polyadics and its reciprocal set 
with respect to multiple dot product are bi-orthogonal. 

The C; may now be eliminated by the aid of (17); and the funda- 
mental identity (15), (on the hypothesis that Co does not vanish), 
rearranged in the form 


M = (A’\A, + A’nAn): M. (19) 


The expression in parentheses will be denoted by J; it is an idemfactor 
in the sense that its A-tuple dot product into an arbitrary K-adic 
leaves that K-adic unaltered; and it is the most important special 
case of a double polyadic,— a dyadic whose antecedents and consequents 
are polyadics. 

Let the expression in parentheses be now transformed by putting 
in the values of all the polyadics A and A’ from the expansions (11) 
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and (18). Collecting the scalar coefficient of E,E; we find it to be 


— Lcjxai, summed on i, which is unity. Again, the coefficient of 


Co 

E,E, when j and k are unequal is summed on 12, which is 
zero. Hence J = LE,E, summed on k, and this latter expression for 
the idemfactor is symmetrical in the sense that it is unaltered by 
interchanging antecedents and consequents. It follows that the 
reciprocal relation between two sets of polyadics is a mutual one, and 
we have 


I = = = (20) 
and also 
M=/:M=M:I (21) 


or the idemfactor may be used either as prefactor or as postfactor. 
To illustrate, if we take K = 1 and N = 3, our polyadics reduce 
to vectors in ordinary space; and if our unit vectors be the usual 
i, j, and k the idemfactor becomes the usual idemfactor of Gibbs, 
namely ii+ jj+ kk. If K = 2 and N = 3 with the same notation, 
our polyadics are ordinary dyadics; the idemfactor (20) is a double 
dyadic, having nine terms obtained by doubling the nine fundamental 
dyads, viz. | 
I = iiii + jjjj + kkkk + ijij + jiji + jkjk + kjkj + kiki + = 
22 
and it may easily be verified that the double dot product of any of the 
nine fundamental dyads ii, ij, etc., either by or into this expression 
gives the dyad unchanged; whence the same is true for any dyadic. 


4. DousLeE POLYADICS. 


Two polyadics A and M written together with no sign between them 
will constitute a double polyad AM. A sum of double polyads is a 
double polyadic, assuming always that the terms of the sum are of like 
character. In the present investigation it will always be assumed that 
the factors A and M of all double polyads are polyadics of the same 
order K. When desirable, all the antecedents A may be expanded in 
terms of the fundamental K-ads E;, and the consequents collected so 
that any double polyadic gy may be written in the form 


y = E,M,+ EM,+---+ E,M,. (23) 
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Each of the consequents M; may be expanded in terms of the E; by 
aid of the notation 


M; = mak, + +- minEn. (24) 


If the M; of (23) be expanded by (24) the complete development of ¢ 
will be in terms of the double polyad units E;E, or 


= 2m (2, k, = 2,° n|. (25) 


It is clear that a double polyadic depends on n?, that is N°“ scalar 
elements. Thus double dyadics in ordinary space, of which the idem- 
factor (22) is a special case, would yield in general 81 terms if it were 
necessary to expand completely, which fortunately it is not. A 
double dyadic is of course a tetradic. If need arose to emphasize the 
2K-adic character of a double polyadic and at the same time to depict 
the scalar elements, use could be made of a system of adjacent squares 
like (6). In the applications which follow, however, the scalar ele- 
ments of a double polyadic are to be thought of as a binary assemblage 
or square array. For example if K = 2 and N = 2, with i and j for 
unit vectors, the order of the fundamental dyads may be agreed upon 
by taking E, = ii, E, = ij, E; = ji and E, = jj. The corresponding 
scalars myq will then be arranged as 


My, M13, Mi4 
Moi, Meo, (26) 
M32, M33, M34 
M41, ‘M42, M43, 


where pq is the coefficient of E,E, in the double dyadic 2m, ,E,E,q. 
Comparing with the four-subscript arrangement (6), we see that the 
upper left hand square of (6) would correspond to the first row of (26), 
the upper right hand square of (6) to the second row of (26), and so on. 

Returning to our fundamental concept of a double polyadic as a 
sum of dyads whose antecedents and consequents are polyadics, we 
next define the A-tuple dot product ¢: ¢g2 of two double polyadics, 
(in precise analogy with Gibbs’ dot product of two dyadics), to be the 
sum of all terms of the form AM : BN or (M : B) AN obtained by mul- 
tiplying out the two double polyadics and taking the K-tuple dot 
product of each consequent of ¢; into each antecedent of go. Thus 
¢1: ¢2 Is also a double K-adic. We abbreviate ¢: ¢ as ¢’, and ¢: 9: ¢ 
as ¢*, etc. This multiplication is easily seen to be associative, as well 
as distributive, but of course not generally commutative. 
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5. Tue Hamitton-Cay.Ley EQUATION FOR DouBLE POLYADICS. 


Let ¢ be a double K-adic expanded as in (23) so that the antecedents 
are the fundamental K-ads E,--E,. Let be a second double K-adic 
> BN and let its consequents first be expanded in terms of E,--E, and 
terms collected so that VY may take the form 


= BE, + BE,+---+B,E, (27) 


so that the consequents are the fundamental K-ads. If we now form 
the K-tuple dot product ¢:¥ it is clear that the scalar coefficient of 
E,E, will be M;:B; or, in terms of the scalar elements of the polyadics, 
will be 2m;.b,, summed on s. This agrees in form with the law of 
multiplication of matrices of order n, a result which might have an- 
ticipated; for if vectors are analogous to polyadics we should expect 
matrices to behave like double polyadics, just as matrices of order NV 
behave like ordinary dyadics. The analogy holds of course only so 
long as we are concerned with formal laws possessed in common by 
the two algorisms. : 

By virtue of this analogy, however, it is evident without further 
proof that any double polyadic ¢ must satisfy an identity of the same 
form as the Hamilton-Cayley equation belonging to a square matrix 
of order n, namely 


mo" + meg" — + (—1)"m,I = 0 (28) 


where the coefficients m,:--m,, are scalars. The determinant of the 
matrix or double K-adic g-gI, where g is a scalar, is 


mi—-9, Me 4, M3, Min | 
My , M3, Man (29) 


and by writing ¢ in place of g and equating to zero we have the Hamil- 
ton-Cayley equation.” 

It is by no means necessary, however, to expand in terms of the 
unit polyads in order to form the Hamilton-Cayley equation. Con- 
sider for a moment the case K = 1, N = 3, when the double polyadic 
becomes Gibbs’ dyadic in ordinary space. The coefficient m, in the 
Hamilton-Cayley equation becomes Gibbs’ ys which, if the dyadic be 
is La-b, or symbolically = a-b. Again, the coefficient mz 
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becomes half (yy) or gos which, if ab and cd are any two terms of the 
dyadic, is the sum of all terms of the form a-be-d — a-de-b, or 
symbolically 


2m. = a-be-d — a-de-b. (30) 


The case of dyadics whose antecedents and consequents are com- 
binatorial products of vectors in N-space, the sum of the classes of 
antecedent and consequent being equal to the dimensionality of the 
space, has also been treated;? and by applying multiple algebra to a 
projective point space of three dimensions, various types of dyadics 
occur, with corresponding Hamilton-Cayley equations.* 

By analogy it is plain we shall expect m, for the double polyadic to 
be the sum of A: M, the A-tuple dot product of each antecedent into 
its own consequent. Using the expansion (23) and noting that E; : M; 
= m; we have TA :M = Ym,; which is the same as m, by the usual 
matrix theory.° If we agree to write, by analogy with Gibbs, SA: M = 
vs, called the scalar of g, we shall have 


m, = ZA :M = gs. (31) 


By analogy with Gibbs’ (¢38)s we may now define the following 
scalar: 


((AM, BN))s=A:MB:N—A:NB:M (32) 


where A, M, B, and N are any four polyadics whatever. If ¢ and @ 
are any two double polyadics the scalar ((¢, @))s is to be found by 
multiplying the double polyadics term by term and adding the scalars 
of all the terms. This scalar will be treated later on as a kind of 
product of the two double polyadics. For the present it will suffice 
to note that if @ = ¢ the scalar becomes twice the coefficient of ¢"~* in 
the Hamilton-Cayley equation, or 


2m. = ((¢, ¢))s. (33) 


The proof is simple. Let ¢ be expanded in the form (23), and note 
that by (24) we have 


E; : My = mui. (34) 


Multiplying ¢ into itself we shall have terms of the form E;M,E,;M, 
whose scalar on developing by (82) vanishes, and terms of the form 
E;M;E,.M;, whose scalar is — tWo-row minor from 
the determinant of the n? elements m,, of g. Each of the latter 
terms will occur twice, with all possible selections of pairs of subscripts 
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Hence the required scalar is twice the sum of two-row minors whose 
main diagonals coincide with the main diagonal of the determinant 
of yg. That is, (by the usual matrix theory), it equals 2me, as we see 
by developing (29). 

More generally, let A;, Ao,---A, and B,, Bo,---B, be any 2p polya- 
dics. Definition. We define the scalar 


((A:B,, A,Bp))s (35) 


to be the sum of terms computed as follows: the leading term is the 
product of K-tuple dot products A;: B,A:: B.---A,:B,; the other 
terms are of like form, and are obtained from the leading term by 
keeping the antecedents A,---A, fixed in position while the conse- 
quents B,---B, are permuted in all ways; the sign of any term is 
positive or negative according as the number of simple interchanges 
needed to form the term is even or odd. 
If g1, g2,°-*, gp are any p double polyadics, the scalar 


((¢1, °°, Yp))s (36) 


is defined to be the result of expanding each ¢ as in (23), multiplying 
out, and adding the scalars formed from each term by the above defini- 
tion. We may now suppose these p double polyadics to be all equal. 
The scalar (36) becomes ((¢, g,--+ to p factors))s and may be abbre- 
viated ((y?))s. It is now easy to see that 


((¢?))s = pimp (37) 


Proof. Expand ¢ as in (23). The terms of the product ¢, ¢,---, to 
p factors are of the form 


where subscripts may have any values from 1 to n; we see by (34) 
that the scalar of this expression has for its leading term mjymj;;mxz 
--+m,,, While, by the definition, the other terms are formed by inter- 
changing the subscripts corresponding to the consequents. Therefore 
the scalar obtained from each expression (38) is a determinant of 
order p whose main diagonal lies on the main diagonal of the determi- 
nant of y, provided the p subscripts 7, j, k, ---, r are all different; if 
they are not all different the scalar of the expression vanishes because 
it is a determinant with two or more identical rows as well as columns. 
Each choice of subscripts of non-vanishing terms will occur p/ times. 
The sum of the scalars of all expressions (38) is then p/ times the sum 
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of p- row minors along the main diagonal of the determinant of the 
miz- By the usual matrix theory, it is therefore p/m >. 

In the above demonstration it is assumed that ¢ has been expanded 
in the form (23). It follows from the distributive character of all the 
steps involved in calculating these scalars that their value is inde- 
pendent of the particular form of ¢, precisely as in the case of similar 
scalars occurring in the usual vector analysis. This is a property of 
fundamental importance. For, when we have to perform operations 
into which these quantities enter, it is frequently sufficient to treat a 
single expression such as (35) in place of the sum of expressions ob- 
tained by expansion of (36); and we may develop ¢, when necessary, 
as YAM in any desired manner. 

To calculate any coefficient m, in the Hamilton-Cayley equation, 
we have then merely to express ¢ as SAM in any convenient way, 
make all possible selections of p different terms, form the scalar (35) 
for each selection, and add the results, dividing by p/. 

Many other analogies with ordinary dyadics will naturally suggest 
themselves. It will be sufficient at this point to note briefly two of 
the more significant of these. 

By solving the Hamilton-Cayley equation as if ¢ were a scalar we 
may find roots 91, g2, -**, gn and corresponding polyadics Ry, Re, ---, 
R,, such that 


— gf): Ri = 0, 1,2,---,n). (39) 


The polyadics R,--R, may by analogy be called the axes of the double 
polyadic ¢ with respect to dot product. 

Again, we may say that a double polyadic which can be expressed 
as the sum of / (and no fewer) properly chosen dyads AM has n—1 
degrees of nullity. To get the simplest criterion of the number of 
degrees of nullity, we may, following Gibbs and Wilson,® introduce 
cross products of polyadics by the law 


AYB = —BYA (40) 


and it is then possible to follow reasoning parallel to that of the paper 
referred to in Note 3, so far as concerns the number of degrees of 
nullity, by use of double powers of ¢. Evidently a product defined 
by (40) is not itself a multiple cross product in the sense that it would 
reduce to double cross product when A = 2. For in Gibbs’ system 
the product by double cross of two dyadics with one common factor 
vanishes, thus abxac = 0, while by (40) the cross product of two 
polyadics vanishes only when they are scalar multiples of one another. 
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It need hardly be emphasized that formal analogies, when they exist, 
are of much value, but do not in themselves imply more than certain 
common laws of operation. 


6. THE ScALaARs (36) ARE SuMs OF CuBIc DETERMINANTS. 


The quantities denoted by (36) are more general than the coefficients 
in the Hamilton-Cayley equation. When p has any value from 2 ton 
such a scalar is a function of all the elements of each of the double 
polyadics which enter. Suppose p = 2 and let ¢; and ~ be expanded 
in the form (23). The scalar elements of g; and yg. may be designated 
respectively by m;; andm:;;. We shall have in the indeterminate 
product ¢i¢g2 pairs of terms of the form E;M,;, E;M.; + E;M,;, 
which, by forming the scalar by the definition (32), gives miyme;; — 
M44 + — These four scalar terms are the same 
as the cubic determinant 

Meoiiy Mei; 
(41) 


M1 ji, M15; Meji, 


where the indices 7 and j are signant,’ but the first index is non-signant 
as might be expected because ((¢1, g2))s = ((¢e, ¢1))g. We see that if 
the respective determinants of gy; and ¢g. be placed as non-signant 
layers of a cubic matrix, the cubic determinant (41) has for its own 
non-signant layers the two-row minors whose main diagonals occupy 
corresponding positions on the main diagonals of g; and g. As a 
verification, if ¢; = g the cubic determinant does not vanish, but 
becomes twice one of these minors, in agreement with the preceding 
discussion of the coefficients in the Hamilton-Cayley equation for ¢. 

Since 7 and 7 have all pairs of values from 1 to n, it is clear that 
((¢gige))g is the sum of all cubic determinants of the second order whose 
main diagonals occupy corresponding positions on the main diagonals 
of the determinants of g; and g. 

An analogous proof holds for all values of p up ton. Thus (36) can 
be written as the sum of all the cubic determinants of order p whose 
main diagonals are corresponding elements from the determinants of 
Ply 

In particular, when p = n, that is, when we have any n double 
polyadics, the scalar (36) is the cubic determinant of order n of which 
the non-signant layers are formed by the square matrices mij, moj; 
*+*, Mnij Of these double polyadics. 
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7. SCALARS FORMED IN COMBINATION WITH THE IDEMFACTOR. 


It may happen that one or more of the double polyadics of the 
scalar (36) is the idemfactor. The scalar formed from any p double 
polyadics together with the idemfactor taken g times is the same as 
the scalar formed from the p double polyadics alone, aside from a 
numerical factor f that is 


((¢1, G2, °° Pps %p))s (42) 


where, however, p + q is not to be greater than n. It is sufficient to 
prove the proposition for one idemfactor and any p double polyadics, 
assuming p+ 1 not greater than n; the general form of (42) then 
follows by letting some of the polyadics equal the idemfactor. I shall 
in fact prove the identity for one idemfactor with p = 2; the proof 
for other values of p is quite similar. The scalar ((¢1, ge, I))g is by 
Art 6 the sum of all cubic determinants of the third order of the form 


1, 0, 0 Miiiy Mik MagK 
0, 1, 0 Mijiy Mizk | Moje (43) 
0, 0, 1 Mikiy Mikj, Mikk Meokj, Mekk 


which on development yields the sum of three cubic determinants of 
the second order having the same form as (41), and corresponding to 
the three pairs of subscripts 27, jk, kt. Consider some particular one 
of these cubic determinants of the second order, say that one which 
involves the subscripts 7j._ It will evidently occur once in the expan- 
sion of every cubic determinant of the form (43), provided the sub- 
scripts are 7, j and any other subscript not greater than n. That is, 
the cubic determinant (41) will occur n — 2 times in the development 
of ((¢1, g2, Z))s. But i and 7 are any pair of different subscripts. 
Hence 


((g1, = (n — 2) ((¢1, (44) 


The same reasoning applies for any value of p. Hence we have 


where p has any value from 1 to n — 1 and ¢:--¢, are any double 
K-adics whatever. By letting several of these K-adics successively 
become equal to the idemfactor we obtain (42) and see that the 
numerical factor f is given by 


p—qt 2)--:(n— p) (45) 


~ 
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We might also have obtained (44) by a count of the number of terms 
obtained by developing either side. If ,C,y,, denote the number of 
combinations of n things p + 1 at a time, this will be the number of 
cubic determinants of order p + 1 obtained on developing the left side 
of (44). Each of these cubic determinants yields p+ 1 cubic de- 
terminants of order p, because one layer is taken from the idemfactor, 
hence contains p + 1 non-vanishing elements. Thus the left side of 
(44) may be written as the sum of (p+ 1),Cp4; cubic determinants 
of order p. On the right the scalar ((¢1, g,-+-, ¢»))s is the sum of 
nC’, cubic determinants of order p. But we have identically 


(p+ (n — p)C> (46) 


agreeing with (44). 

In particular it follows from (42) that any scalar (36) can be written 
as a single cubic determinant of order n by adjoining the idemfactor 
n — p times, so that p+ q = n and the numerical factor becomes 
(n — p)!. Thus 


((¢1, $2°°*s Pps (n p)! ((¢1, Yp))s (47) 


in the case where the idemfactor occurs just n — p times. 

It also follows that any coefficient in the Hamilton-Cayley equation 
for g can be written as a single cubic determinant of order n. For by 
letting the p arbitrary double polyadics all equal yg, we have in virtue 
of (37) 

I, I,---1))g = p! (n — p)! my (48) 


provided ¢ enters p times and J enters n — p times. 


8. ScaLARS EXPRESSED AS SQUARE DETERMINANTS WHOSE ELE- 
MENTS ARE K-ADICs. 


By virtue of the result of Art. 7, it becomes possible to express any 
scalar of the form (36) without the use of cubic determinants, if we 
are willing instead to use determinants whose elements are not scalars 
but polyadics, and to make certain slight but needful modifications in 
the usual rules for determinate expansions. 

Suppose, as before, that there are n arbitrary double polyadics 
¢1°**¢n and let each of them be expanded in the manner of (27), thus 
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where the consequents E,---E, are the fundamental K-ads and the 
antecedents are any arbitrary K-adics. Let E; be any one of the 
fundamental K-ads. We have 


gi: EB; = (50) 


The binary assemblage of n? arbitrary polyadics B;; are now to be 
formed into a determinant 


Bu, By, Bin 
Boi, Boo, Bon (51) 


B..2, Ban 


where the law of expansion is defined to be as follows: the leading 
term is written [B,,B.---B,,| and signifies the determinant in the 
ordinary sense whose elements are the scalar elements of the polyadics 
which enter: namely, if we expand as 


By; = bi + +++ ++ (52) 


the leading term in expanding (51) is the determinant of ordinary 
character 


bin, bir, Din 


boo0, boon (53) 


Dants Dann 


the other terms in the expansion of (51) are of similar character, and 
are formed by writing out the determinant (51) as if the elements were 
scalars, enclosing each term in square brackets to signify the determi- 
nant formed from it as in (53). The signs will occur as in ordinary 
determinants if we adopt the rule that, in every term, the order of the 
polyadics shall be the order of the rows in the original determinant. The 
meaning of the determinant (51) is thus defined without ambiguity. 
We note that by this rule the order of rows in (53) corresponds to the 
order of rows in (51), in the sense that each row of (53) consists of the 
elements of a single polyadic from the corresponding row of (51). 

It is then apparent that (51) is the sum of n/ ordinary determinants, 
hence the sum of (n/)? terms on complete expansion. The order of the 
rows is non-signant, in the sense that, if two rows in (51) be inter- 
changed, the value of the expression is quite unchanged; for on ex- 
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panding, first, into a sum of ordinary determinants or square bracket 
expressions, the sign before every square bracket is changed because 
the signs behave as if (51) were an ordinary determinant; but, since 
the development follows the rows, the sign of every determinant (53) 
will have been changed due to the interchange of rows; thus each 
square bracket has its sign changed, we may say, both from within and 
from without. Hence the whole expression is unaltered. On the 
other hand, an interchange of two columns in (51) changes the sign 
before every square bracket but leaves the brackets themselves un- 
affected, for the development is according to rows; hence the sign of 
the whole is changed. 

It is now not hard to see that the polyadic determinant (51) is the 
same as ((¢1, ¢2,°**,¢n))g. For we have seen that the first index is 
non-signant, the second signant; the third is obviously signant since 
it governs the columns of the ordinary determinants of the form (53) 
but does not otherwise enter. Furthermore, on complete expansion 
of (51) we shall have each term a product of n elements, one from each 
of the double polyadics, one from each row of (51), and one from each 
column of (51). But this is precisely the law of formation of a cubic 
determinant having the elements of the polyadics as non-signant 
layers. By Art. 6 this is the scalar in question and we may write the 
identity 


g:E, gi:tE, 
¢g2: Ey, 2: En, 


((¢1, Yn))s (54) 


where the polyadic determinant on the left is defined as above. 

To take the simplest possible illustration, let K = 1, N = 2,n = 2. 
Our K-adics become vectors in two dimensions. Take E, = i, E, = j. 
Let there be two dyadics gi = + Dioj = biti + ij 
by20jj and similarly for g. The determinant on the left of (54), be- 
comes 


whose elements are vectors. Expanding by the rule we have 
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which is the same as the difference of two determinants 
(Dy11b222 — bi12b201) — — bi22be11) 


and this, again, agrees with the cubic determinant formed with the 
element of the two dyadics as layers, namely 


beri, 
bore, boo» 


A number of results follow immediately. We might have formed 
the left side of (54) by taking the multiple dot product of E,, Es, ete. 
into each of the arbitrary double polyadics, as E,: ¢g;, Ez: ¢, etc. in the 
first row, and similarly for the other rows. Instead of (51) we should 
have had the determinant of the M;; by developing after the manner 
(23) instead of (27). The rest of our reasoning would have been un- 
changed, however, hence the value of the expression would have re- 
mained the same. The change would in fact be equivalent to chang- 
ing rows into columns and columns into rows in all the non-signant 
layers of the cubic determinant, which, since the second and third 
indices are signant, leaves the result unchanged. This is illustrated 
by the example just given, where the last index refers to the rows in- 
stead of to the columns as in (41). If we let y’ stand for SMA when 
yg stands for SAM, and say that the conjugate ¢’ of a double polyadic 
y is formed from it by interchanging antecedent and consequent, it is 
clear that E:g = ¢’ :E. We therefore have the identity 


bin, 


which may also be deduced directly from the definition of these scalars. 
The well known fact that the Hamilton-Cayley equation for any 
dyadic or matrix is identical with that of its conjugate appears as a 
special case of this identity. 

Also if n — pout of the arbitrary double polyadics are allowed to 
become equal to the idemfactor, the identity (54) becomes in virtue 
of (47), 


gi: gi: Ep, ’ gi: E,, 
go: go: Ey, 2: E, 
¢,: By, Ey, E,, 
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where the last n — p rows are alike, a circumstance which need cause 
no surprise, since, as already pointed out, the rows are non-signant in 
the very nature of the definition of these polyadic determinants.® 

If we still further specialize (56) by letting the remaining p double 
polyadics all become equal to ¢ we shall have, by (48), 


¢ :E,, E,, E, 

E,| = p! (n — p)! Mp (57) 
E,, 


Ei, E, 


where, on the left, p rows are like the first row, and (n — p) rows are 
like the last row. 


9. AN INVARIANT PROPERTY OF THE SCALARS (36). 


From the form of the definition of the scalars of the present discus- 
sion, they are independent of development in terms of a particular 
set of polyadics. Throughout the argument the fundamental polyads 
E,--E, may be replaced by any other set of » polyadics forming a 
normal, orthogonal system with respect to K-tuple dot product, for 
no other properties of these polyads have been so far used. These 
scalars have another kind of invariance of a somewhat wider sort, with 
respect to a reference system which need no longer be normal nor 
orthogonal but only linearly independent. The properties in question 
are essentially included in the following theorem: — 

Theorem. The scalar ((¢1, ¢,--++,¢n))s formed from a set of n 
arbitrary double polyadics is unaltered when each of these double 
polyadics is multiplied by (or into) the same double polyadic 6, (by 
K-tuple dot product), except for a factor which is the determinant of 6. 

This theorem is equivalent to the identity 


G2, $n))s ((6"))s = n! ((¢1: +++, On: 9))s (58) 


for, by (37), ((0"))s equals n/ times the coefficient of 6” in the Hamilton- 
Cayley equation for 0, which is the determinant of 6; and 6: g = 
gi: 9, while by (55) the left side of (58) is unchanged if we change 
every ¢ into yg’; hence it is sufficient to prove (58) as written. 
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Proof. Let every ¢ be expanded as in (23). The scalar ((g: --, 
¢n))s is equal to the polyadic determinant of the M;; as shown in 
Art. 8. To multiply every ¢ into a double polyadic @ is equivalent to 
performing the same linear transformation upon every M. But the 
polyadic determinant is the sum of ordinary determinants whose 
rows consist of elements of the respective polyadics M. It is well 
known that when every row of a determinant is affected by the same 
linear transformation the determinant is merely multiplied by the 
determinant of the transformation. Thus every square bracket 
term in the expansion of the polyadic determinant of the M;; is multi- 
plied by the determinant of @. The theorem is therefore proved. 

As an immediate consequence, any polyadic element ¢;:E; on the 
left of the identity (54) may be replaced by ¢;:6:E;, provided we 
multiply the right side by the determinant of 0, which may be called 
mg and is the same as the determinant of the n polyadics 0: Ej, 6: EB, 
--,0:E,. By a proper choice of the arbitrary double polyadic 8, 
these n polyadics may be made to be any n polyadics whatever, Aj,---, 
A,. If we agree to write [y;:A,] for the polyadic determinant whose 
elements have the form indicated, and [a;;] for the determinant (in 
the ordinary sense), whose elements are the scalar elements of the 
polyadics A,--A, resolved according to (11), we shall now have the 


identity 


[ei: Ai] 
((¢1, $n))s = 69) 


assuming the polyadics A; linearly independent. Thus the — on 
the left is invariant of the particular choice of these polyadics.® 
In a similar manner and with similar notation we may derive from 


(56) 


Aj, go: A;, Pp: A;, A;, A;, , Aj] [a;;] (n—p)! 
((¢1, (60) 


and from (57) 
A;)? (A;)"-?] = [ai,] p! (n — p)! mp (61) 


In employing this abbreviated notation for polyadic determinants, if 
there are two subscripts the first always refers to the rows, the second 
to the columns, as in (59) and (60). Where a series of n polyadies is 
written in brackets as in (60) and (61) these refer to the rows. We 
pick out an element of a particular column by giving a particular value 
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to the common subscript 7. If several rows are alike the fact may be 
denoted by an exponent, as in (61). Thus in (60) the left side might 
have been written (A;)"~?]. Where, however, the 
notation implies merely n polyadics (instead of n”) the square bracket 
signifies the ordinary determinant of the n? elements by the scheme of 
(11). Thus instead of [a;;] we might write [A;] without change of 
meaning. Unless the contrary is stated, all subscripts are to run 
from 1 to n. 


10. Forms oF THE SCALARS (36) WHICH SHOW THEIR PoLyapic 
CHARACTER. 


Up to this point have been noted three distinct methods of actually 
calculating a scalar of the form (36), namely 

1°. Multiply the double polyadics term by term and take scalars 
according to definition. 

2°. Form a sum of cubic determinants. 

3°. Form a polyadic determinant by operating on any set of 
linearly independent polyadics and divide by the determinant of the 
set. 

In all these cases the double polyadics have behaved formally as if 
their antecedents and consequents were vectors in space of n dimen- 
sions, instead of polyadics each of which is a sum of products of K 
vectors in space of N dimensions, where n = N*. Each of these 
methods applies with no change whatever if K = 1 so that n = N and 
the polyadics become actually vectors,— with the obvious exception 
that the multiple dot product, indicated by a colon, becomes ordinary 
or single dot product. While I am not aware that any of these three 
methods has actually been employed by writers on vector analysis in 
N dimensions, yet it can hardly be said that the individual steps do 
not occur, in principle at least, in the writings of Hamilton, Gibbs, or 
Grassmann. For example, the definition of these scalars is equivalent 
to a double indeterminate product followed by inner product. In the 
case K = 1, N = 3 they were introduced by Gibbs as scalars of double 
cross products. In the same case they were expressed by Hamilton 
both by summations which are equivalent to double multiplication,’ 
(though never so regarded by him), and also as quotients equivalent 
to our method 3°. It is true that none of these writers expressed the 
scalars as sums of determinants whose elements are extensive magni- 
tudes, use of which was first made by Cayley, nor by cubic determin- 
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ants, also due to Cayley, but these two equivalent formulations differ 
more in appearance than in reality from such expressions as Hamilton’s 


= Sigigk + Sipkgi + Skyigj 
which in Gibbs’ notation would be 
Me = 


where 2, j, k denote rectangular unit vectors. 

It will be of some value for the sake of later applications if we take 
account of these scalars as they appear when our polyadics are de- 
termined by K-dimensional matrices. Consider for a moment the 
case ((AM, BN))s where A, M, B, and N are dyadics in two dimen- 
sions,— evidently the simplest case which can be regarded as a gener- 
alization of the usual vector analysis. With unit vectors i and j 
we may either write 


A = aylii + apvij + anji + anjj (62) 
or may regard A as an abbreviation for the matrix 


ay, Ay 


Ar (63) 


with similar notation for the other three dyadics. If we form the 
required scalar directly from the definition, namely as A: MB: N — 
A: NB: M we obtain 


((AM, BM))s = (aumu + + + (bunn + 
+ boyne: + 
+ boome2) (64) 


If instead we use the method of cubic determinants we shall first 
form the square array of the fourth order belonging to the double 
dyadic AM, which is 


a2) 
11, 


@22™21, 


A222 


The first of these will be 


(65) 


with a similar array for BN. We then take the sum of all cubic 
determinants along the main diagonal, using the two arrays as simi- 
larly placed layers. 
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bunu, bine ( 6 6) 
which develops into 
— — + (67) 


There are six of these groups of terms, or 24 in all. The scheme of 
each group is given by 

Apgm | pqit (68) 


ArsMpq, 


which develops into | 
— ArsMp gNrs — ArsMrDpgNpq (69) 


and the summation is performed by giving to the number-pairs pq 
and rs unlike values chosen from the pairs 11, 12, 21, and 22. Com- 
paring with (64) it is easy to verify that the two methods agree. 

If we use the method of polyadic determinants and adopt the dyads 
ii, ij, etc. as operands we first form the dyadics AM: ii, AM: ij, etc. 
treating BN in a similar manner. We thus have two rows of a de- 
terminant of the fourth order whose elements are dyadics. The other 
two rows are alike and given by /: ii, J: ij, etc. that is by ii, ij, ete. 
The result 


Amy, Amz, Amn, Amz 
Bra, 
ii ij 


should by (56) be the double of the required scalar. Developing by 

the definition of Art. 8 the leading term is [Amy, Bny, ji, jj] which 

denotes the ordinary determinant of the fourth order 

AyMy, AyMy, 

(71) 


which develops into aymybyn~2— aymybyn». The next term may be 
taken as — [Am, Bny, jj, ji] which doubles the two scalar terms 
already found. It is easy to see that developing (70) by two-row 
mjnors after Laplace’s method yields the scalar terms in a fashion 
similar to the method of cubic determinants. 
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The fact that any scalar (36) is unaltered by changing all the double 
dyadics into their conjugates, and is also unaltered by changing the 
order of double dyadics, appears in the double symmetry of (64). A 
further and quite distinct fact could not have been brought out so 
long as our notation represented polyadics as formal vectors: a scalar 
formed from double dyadics, like (64), is unaltered if all the dyadics 
A, M, B, N are replaced by their conjugates. This is equivalent to 
interchanging the two subscripts throughout the expression on the 
right of (64); whereas changing the double polyadics into their con- 
jugates is equivalent to everywhere interchanging a with m and at the 
same time b with n; and changing the order of the polyadics is equiva- 
lent to everywhere interchanging a with 6 and at the same time m 
with n. 

More generally a polyadiec which is a sum of terms of the form 
@&-*:a&, where the K factors are vectors in space of N dimensions 
may be regarded as one of a set of K! polyadics obtained from one 
another by making like permutations of factors in every term. A 
glance at the definition (10) of multiple dot product is sufficient to 
show that all scalars (36) are unaltered when all the antecedents and 
consequents of every double polyadic have the factors of all their 
terms permuted in the same manner. If each polyadic consists of a 
single polyad and we apply the definition of Art. 5, the truth of the 
proposition is evident, hence by the distributive principle holds 
universally. 

Examples like (64) or (69) may be generalized in three ways, ac- 
cording as we increase N, the dimensions of the space, K, the order of 
the polyadics, or p, the order of the scalars. 

If N increases we still have double dyadics and the form of (69) is 
unaltered, but the summation is performed over a larger group of 

number-pairs, namely 11, 12, 21, 13, 31,---, NN. 

If K increases the number of subscripts i increases but (69) still con- 
sists of 4 scalar terms and the summation is over pairs of unlike 
number-triplets, quadruplets, etc. 

If p increases the order of the determinant (68) increases, likewise 
the number of terms in (69). Thus if we keep N = 2, K = 2, but 
make p = 3 there are 4 cubic determinants whose sum is ((AL, BM, 
CN))s or 72 scalar terms in all. 

If we have double dyadics in 3 dimensions with p = 2, ((AM, BN))s 
will be the sum of 36 expressions of the same form as (69) or 144 scalar 
terms. 

It is usually possible to neglect the dimensionality of the space and 
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the number of terms in the double polyadic, frequently also the order 
of the polyadics. Thus ((AM, BN))s is symbolic of ((¢1, ¢e))g. 
When the order of the polyadic is given we may also treat with polyads 
in many operations as symbolic of polyadics. Thus if we have double 
dyadics in any space we may often write aa’ for A, mm’ for M, etc. 
in which case a bar may be used to separate antecedent from conse- 
quent in the double polyads, or ((aa’| mm’, bb’| nn’))s may be taken 
symbolically for the much more general ((¢:, ¢2))s. Hence if g and 
¢2 are double dyadics we may write 


((¢1, ¢2))s = 2(a-ma’-m’b-nb’-n’— a-na’-n’b-mb’-m’) (72) 


and may even omit the summation sign in many operations. It is of 
interest to compare this vector expression with (64) which was based 
on the notion of the dyadic as a matrix. 


11. INVARIANTS REGARDED AS PRODUCTs. 


The examples of the last article are sufficient to indicate that the 
scalars (36) are invariants of polyadic systems and that they are 
subject to a considerable variety of transformations. There is one 
more method by which they may be calculated, based on a concept of 
these scalars which may serve to widen considerably our view of the 
formal laws to which double polyadics are amenable. Returning to 
the earlier conception of a double polyadic as a sum of unit double 
polyads E;E; each multiplied by a scalar, as in (25), consider scalars 
of the form ((E,E;, E,E,))s, that is E,;:E,E,:E,— E,;: E,E,: Ex. 
The unit polyads E; and E, are antecedents while E; and E, are conse- 
quents. We see that the scalar is equal to +1 when each antecedent 
is equal to its own consequent, but unlike the other antecedent and 
consequent; is equal to —1 when each antecedent differs from its own 
consequent but is equal to the other consequent; and is otherwise 
zero. That is 


((E,E;, E,E.))s= +1, ((E:Ex, E.E;))s = — 1, ((E,E;, E,E,))g = 0 
((E.E;, E.E;));=0, ((E,E;, = 0 
E.E,))s 0, ((E;Ex, E.E,))s = 0 


(73) 


on the understanding that subscripts are unequal unless denoted by 
the same letter. 
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In the same way it is easy to see that any scalar derived from a set 
of p unit double polyads, namely 


((E,E,, E,E;, E,E.))s (74) 


where subscripts run from 1 to n, will vanish unless, by interchanges 
among the consequents, it can be brought to the form 


((E,E,, E,E,, E,E,))s (75) 


wherein each antecedent is equal to its own consequent; and moreover 
vanishes unless all antecedents are unequal. The sign is plus or minus 
according as the number of interchanges needed to give the form (75) 
is even or odd. 

We may now regard the scalars (36) as products of double polya- 
dics formed according to this law, and, if each double polyadic is 
developed as in (25), can evidently calculate these scalars with no 
more labor than would be needed to multiply algebraic polynomials. 


12. Srar Propvwcts. 


It is especially in the case p = 2 that the concept of the last article 
appears useful. In keeping with the idea that ((¢1, ¢2))s, the calcula- 
tion of which has been exemplified in some detail, is a kind of product 
of the two double polyadics, we may introduce the notation 


((¢1, ¢2))s = 91" (76) 


and speak of this scalar as the star product of the two double polyadics. 
It is essentially a kind of double product, but since there is nothing 
analogous for mere vector factors it is not necessary to write two stars 
to emphasize the distinction. It may properly be compared, or rather 
contrasted, with the double dot product of Gibbs. To illustrate, take 
again the simplest case, K = 1, N = 2, so that the factors are dyadics 
in space of two dimensions. By Art. 11 the multiplication table for 
unit dyads will be 


ii, ij, ji, jj 
@ @ +41 


& Gai, (77) 
@& 
jj + k 0, 0, 0 
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which brings out the ‘‘skew” character of the star product: it is 
formed by the multiplication of unlike elements. If the scalar ele- 
ments of ¢; and g are a,, and b;, in the style of (63), we shall have 


(ayii + + anjj)*(buii + disij + doiji + beojj) 
= — — + (78) 


The result is of course the same as if the dyadics had been expanded in 
any other manner and we had applied one of the three methods of 
Art. 10. 

By contrast the double dot product ¢:: ¢2 yields a multiplication 
table having +1 along the main diagonal and elsewhere zero, giving 


$1: = Anbu + + + (79) 


The concept of ¢:*¢. as a product is justified by the fact that it is 
in every case possible to set up systems of double polyadics ¢, 
°*, such that 1 if subscripts are equal, otherwise zero. 
Here v = n? = N* for this is the number of linearly independent 
double polyadics in any case. In the example (77) if we make 


(ii + jj) (ii — jj) (ij + ji) (ij — ji) 
’ ’ F;= — F,= 


gi" 


F,= 


V—2 


we shall have F,*F,= 1 if subscripts are equal, otherwise zero. 

It is clear that, just as the dyads ii, ij, ji, jj could be formally 
treated with reference to double dot product as would be vectors in 
4-space with reference to ordinary dot product, so the dyadics F can 
be similarly treated with respect to star product. 

A correspondence may be thus set up between invariant properties 
of the dyadics or polyadics or double polyadics and the geometry of 
the space in which these magnitudes behave formally like vectors. 

A special case of much interest has been treated in great detail, with 
some differences of terminology, by E. Waelsch.44 By his very in- 
genious system the whole of the well-worked theory of invariants of 
binary forms is employed as algorism for phenomena of various sorts 
in ordinary space. In the language of the present discussion his 
system would arise by making up three fundamental units 


+i) 
(81) 


so that E,, E,, E; form a normal orthogonal system with respect to 
double dot product. A binary quadric may be written 


E, = ii, E, = 
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yay? + 2 + (82) 
and a variable vector may be written 
x= Xi Xoj (83) 


but this vector is regarded as merely formal. We have identically 
[ay ii + a2(ij + ji) + EX = + + (84) 


so that any binary quadric is the double dot product of a dyadic into 
the dyad xx. We now make up the new system, orthogonal with 
respect to star product. 


and also write 
_ O11 _ ~ 


when we shall have identically 
+ ar (ij + ji) + (87) 


The system Fj, F:, F; is now taken as a set of rectangular unit vectors 
in ordinary space. Thus any binary quadric corresponds to a vector 
drawn from a stipulated origin in space. As a single but important 
example of very numerous correspondences which arise, we have 


y+ 2= 2 (anax — ay’) (88) 


This brief outline renders scant justice to the very beautiful system of 
Waelsch, from which it differs in slight details, owing to difference in 
definitions and point of view. It differs in particular by identification 
of the dyadic with the desired vector through the star product and the 
relation (87). The essential thing is the correspondence set up be- 
tween geometry and invariant theory by relations of the type (88). 

More generally, if we take the ordinary dyadics of Gibbs, making 
K = 1, N = 8, with unit vectors i, j, k, it is easy to verify that the 
system of nine dyadics given by 


F.V2 =ii+jj, FeV —4=kk—2ii, —4=2jj—kk, F,V—2= 
if ij+ ji 

F;V —2=jk+kj, Fy —2=ki+ik, FiV2=ij—ji, FsV2=jk—kj, 
F,V2=ki—ik 
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are normal and orthogonal under star product. As already noted, 
the star product in this special case of two dyadics ¢ and @ is the scalar 
(y30)s of Gibbs. If we wish to restrict ourselves to self-conjugate 
dyadics we have merely to drop Fy, Fs, and Fy, from the list. It would 
be interesting to work out the analogy between the scalars called 
above ((¢?))s, which might be symbolically written ((ab))?s, and the 
invariants of ternary quadrics written in the usual symbolic language. 
We see however that systems orthogonal under star product will 
always be possible, because to every double polyad E,E, corresponds 
always at least one other double polyad not orthogonal to the first. 


13. SCALARS FORMED BY STAR MULTIPLICATION. 


In what follows I shall, for simplicity of language, take K = 1, so 
that double polyadics become dyadics, and polyadics become vectors. 
This imposes no formal limitation, and can be extended to the most 
general case with slight verbal change. 

Supposing N to have any value, we assume F,, F.,---F, a set of 
linearly independent dyadics orthogonal under star product. Any 
other dyadics A, M, B, or N can be expressed in terms of these. 
Let ¢; and ¢g2 as before be a pair of double dyadics SAM and 2BN. 
The star product of AM into BN is defined to be A(M*B)N which is 
the same as ANM*B. Thus ¢,*g is a new double dyadic. If we 
regard y*y and g*y*y etc. as a new type of power abbreviated ¢*, 
y**, etc. no further proof is required to show that ¢* satisfies a Hamil- 
ton-Cayley equation similar in form to (28), 


my* (1) (—1) m*, I" = 0 (89) 


nor to show that the star idemfactor J” is identical with =F,;F,;. The 
whole of Art. 3 might be repeated at this point, replacing multiple dot 
product by star product and E by F. The coefficients m*, in (89) 
will be formed as in Art. 5, replacing multiple dot by star product. 
We shall also have a set of scalars 


(¢1, (90) 
analogous with (36). For example 
(AM, BN).*= A*MB*N — A*NB*M (91) 


which differs from (32) in having star in place of colon. It also differs 
in having no analogue when A, M, B, N are mere vectors, whereas 
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(32) would still hold if single dot be written for colon. As before re- 
marked, the star product is essentially double. 
For further comparison, suppose NV = 2, with the notation of Art. 10. 


(AM, BN)*s — — 11) bo Ni2 


+ beens1) 
— — — + (Dyym22 — 
— boymy + beom1) (92) 


which, like (64), reduces to 24 scalar terms by expanding. 

The method of cubic determinants might be applied to star scalars, 
but with less ease, because a system of dyadics orthogonal under star 
product is in general imaginary. 

Identities like (44), (42), (45), (47), (48), and (54) to (61) inclusive 
may be proved in quite parallel steps, with appropriate change of 
notation. With exception of (54), (56), and (57), these identities do 
not imply the presence of V — 1 even when we change to star product. 
While the theory of invariants lies outside the scope of this paper, it 
may be noted that star scalars, like those of (36), are by nature of their 
definition, invariant under transformation of the fundamental units, 
and some of them are invariant in a much wider sense. 


14. IDENTITIES WITH THE STAR IDEMFACTOR. 


If we have a system of linearly independent dyadics A, As,---A, 
their reciprocals with respect to star product may be called A’, 
A”,,---,A”,. By steps as in Art. (3) we may show that DAA” is 
identically the idemfactor 2FF. It is evident by Art. 12 that the star 
idemfactor is real. Taking N = 3 with unit vectors i, j, k, it is easy 


to verify that the star idemfactor 


1 
"= 5 ii(jj + kk — ii) — (ijji+ i) (93) 
because I’’* ii = ii, J’’* ij = ij, etc. In fact since, as operators, 
['*= |: (94) 
both being identical operations, we must have 
A” .*= A’,: (95) 


where A”. and A’, are respectively the star and dot reciprocals of A, 
in the set of linearly independent dyadics Aj, Ao,---,A,. Operating 
with both sides of (95) on the dot idemfactor gives 
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A’. (96) 
or in words: the dot reciprocal is the star product of star reciprocal 
and dot idemfactor. Similarly by operating on the star idemfactor, 

A” (97) 
or: the star reciprocal is the dot product of dot reciprocal and star 
idemfactor. Thus the star reciprocal of ii is ii: J’ or } (jj + kk—ii), 
and the star reciprocal of ij is —ji etc. 

Star multiplication may always be replaced by dot multiplication 
by introducing the double dyadic /*7 between the factors. For if M 
and N are any two dyadics we have M: J = M and J: N = N, hence 

M*N = M:/*/:N (98) 
We also observe that /*/ is the reciprocal (double dyadic) of I’”’ with 
respect to double dot multiplication, that is 

a = I (99) 


In a similar manner double dot multiplication can be changed to star 
by the operation J’: I” for 


M:N = :/'*N (100) 
and I’’: J’’ is the reciprocal of J with respect to star multiplication, 


It is thus plain that, in spite of the imaginary character of star 
orthogonality, the reciprocals of a set of real dyadics are real, and the 
reciprocal, (or inverse), of a real double dyadic, if it exist, is also real. 


15. ILLusTRATION oF Various Types OF MULTIPLICATION. 


We may close this discussion of the general properties of double 
dyadics by a few examples of the types of product considered. Con- 
sider the two tetrads ijji and jiij. If they be taken as polyadics of 
order 4 and multiplied by multiple dot product we have 


(ijji) : (jiij) = 0 (102) 


where colon means multiple dot product. 
If they be taken as a pair of double dyadics the scalar (36) is 


((aj | ji, Ji] ij) = ij: jijis ij — ij: ijji: ji = —1 (103) 


where the colon stands for double dot product, for we assume the case 
K = 2, hence might also have written (ij | ji)*(ji| ij) for this scalar. 
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If we assume K = 1 we return to the case of Art. 13 and 14 and have 
(ij | ji, ji | ij)*s = ij*jiji*ij — ij*ijji*ji = +1 (104) 

and eiisy if we had merely four dyads ij, ji, ji, and ij, the scalar 
((ij, ji, ji, ij))s = 0 (105) 


since all such scalars vanish unless the antecedents are distinct. 

Again, we may form double and mixed products of double dyadics. 
Assuming K = 1 the star product written singly will refer to the 
dyadic factors. A double star will refer to antecedents and conse- 
quents as in the usual double product, thus if A, M, B, N are dyadics 


AM**BN = A*BM*N (106) 
and we may have 

AM*: BN = A*BM:N (107) 
and 

AM:*BN = A: BM*N (108) 
while 

AM:: BN= A: BM: N (109) 


It is evident that, regarding the double dyadic as an element, we 
could start a fresh investigation of each of these types of product with 
respect to orthogonality, idemfactors, and so on ad inf. But we have 
now ample tools for the applications it is desired to make. 


PART II. THE SOLUTION OF THE LINEAR MATRIX 
EQUATION. 


16. TRANSFORMATION FROM MATRICES TO DouBLE Dyapics. 


Any matrix A of order N is defined as a binary assemblage of ele- 
ments a;, where subscripts run from 1 to N. We add two matrices 
by adding corresponding elements. We multiply by the rule 


(AB) x= sx (s 2,° N) (110) 


If Aj, Ao,-++, A, and B,, Bo---,B, and C are known matrices 
while x is a required matrix, the linear matrix equation may be written 


O(a) = AyrB, + = C. (111) 
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The number of terms hf on the left is said to denote the eztent of the 
operation @ or of the equation. 

It is known that any matrix A is equivalent to a dyadic A or Daa’ 
where a and a’ are vectors. Any term A zB is therefore equivalent 
to a sum of terms of the form aa’ -X-bb’ where a, a’, b, and b’ are 
vectors in space of N dimensions and X is a required dyadic. Since 
by the definition of double dot product we have ab: xy = a-xb-y, 
and since X might be written 2xy we shall have 


aa’ X bb’= ab’a’b: (112) 


which is the transformation at the basis of the method to be exhibited. 
Thus the linear matrix equation is equivalent to 


g: X= =MN:X=C (113) 


where ¢ is a dyadic whose antecedents and consequents are themselves 
dyadics, that is a double dyadic. 

Since g: X = 6(x), 6 and ¢ satisfy the same Hamilton-Cayley equa- 
tion. One method of solving the linear matrix equation is therefore 
by putting 


— — +- ++ (—1)"1C (114) 


where m,:--m, are the same coefficients as in (28) and where mz is 
assumed not zero. 


17. TRANSFORMATION BACK TO MATRICES. 


If the matrices A and B are known in such form that transformation 
to (113) is convenient, any of the methods exhibited in Part I may be 
used to calculate the coefficients m. <A general solution, however, 
demands ability to compute them directly from A;--A, and B,:- By 
without the necessity of first forming ¢ by the rule (112). In other 
words we need to transform the scalars (36) from the language of 
double dyadics to that of matrices, taking account of (112). This 
may be conveniently accomplished by a partly symbolic notation, 
omitting subscripts and summation signs which refer to the terms 
1, 2,---A of the matrix equation but occasionally preserving those 
which refer to the p/ terms of the expansion of (35). We thus write 
symbolically 


y = ab’|a’b = MN (115) 
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and m, = ab’: a’b 
a: b’ 
AsBg (116) 


where Ag and Bg are respectively the first coefficients in the Hamilton- 
Cayley equations for A and B; or, what is the same thing, they are 
the sums of elements along the main diagonal of these matrices. 
Stated in complete form the result is 


m= AigBig + +-+-+ AngBras (117) 


As remarked earlier, the symbolic process is valid by virtue of the 
distributive character of all the steps involved. Thus for m: we may 
put, in abbreviated language, 


2m2= ¢))s = ((ab’ | a’b, ed’ | c’d))s 
= ab’: a’bed’: c’d — ab’: c’ded’: a’b 
= a-c’b’-dc-a'd’-b (118) 


Now the first term on the right is the product of the scalars of the four 
dyadics or matrices which enter; that is, it is symbolically AgBsCsDg 
or AigBigAjgBjs. The second term may be further transformed 


a-c’b’-dc-a’d’-b = [a-c’a’-c] [b-d’b’- d] 

[aa’-ce’|s [bb’- dd’|> 

[AC]s [BD]s 

(A :A [BB ils (119) 


where [A ;A,]g is the first coefficient in the Hamilton-Cayley equation 
for the matrix A ;A , that is, for the product of A; and A,, it is the sum 
of elements on the main diagonal of this product; and similarly for 
[B:B;]g. Stated in complete form the result is 


2m: = isBisAjsBjs — [t, 7 = 1, 2,--, A] (120) 


the summation being now with respect to all possible pairs of terms 
selected from the left side of the original equation (111). It is to be 
particularly noted that the summation includes the cases 1 = 7; in 
other words, by “ pairs of terms” is meant not merely different terms 
but also the same term taken twice. The reason is not far to seek. 
By the transformation (112) any term AxB was changed into a sum of 
terms of the form MN: X, but not, in general, into a single such term. 
It follows that, although ((MN, MN))s, which is the star product of a 
single term MN into itself, must vanish, yet a scalar A*g B’s — 
(A?)9(B?)s symbolically derived from it does not in general vanish. 
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This fact will be fully exemplified below. A symbolic equation 
implies, not the equality of individual terms each to each, but equality 
of the sums of all possible terms of the equated forms. 

Following a similar procedure to find m3 we have 


6m3= ((¢, 9, ¢))s= ((ab’ | a’b, cd’ | c’d, ef’ | e’f))s, symbolically, 
a-a’'b-b’[c-e’c’ -e| [d’-fd-f’] 

— bf-b’] — e-e’f-f’[a-c’a’-c] [b’-db-d’] . 

+- (ab’: c’d) (cd’: e’f) (ef’: a’b)+(ab’: e’f) (cd’: a’b) (ef’:c’d) 

(121) 

by direct application of the definition of Art. 5, and by grouping the 
dot products in the second, third, and fourth terms after the manner 
of (119). The first or leading term is evidently the same symbolically 
as AigBigAjgBjgArgBis. The next three terms taken together are 
symbolically — 2A igBig(A;Ax)s(B;Bx)s. In transforming the fifth 
term the procedure, as in all cases, will be to group into one factor 
those vectors which correspond to prefactors A in the original equa- 
tion, and into another factor those which correspond to the post- 
factors B. Thus 


(ab’: c’d) (cd’: e’f) (ef’: a’b) = (a-c’c-e’e-a’) (b’- dd’ -ff’-b) 
= (aa’-ee’-cc’)s(bb’- dd’ - ff’), 
(A (122) 


It is especially worthy of remark that the order of the matrices in one 
factor is the reverse of that in the other, a consequence of the trans- 
formation (112). When there are only two matrices in each factor, 
as in (119), it was of no importance whether we wrote (AC)g or (C/A)g, 
for these are equal: the scalar of the product of two matrices is inde- 
pendent of their order, because (aa’-cc’)s = a-c’a’-c = (cc’-aa’). 
But by similar reasoning we see that the scalar of the product of several 
matrices depends on their cyclic order, as is well known. 

In the same way the last term of (121) becomes (A;A ;Ax)5(B;B.B;)s. 
Collecting results, the complete statement for ms is 


6m3 = A jgAxgB igB jgBrg — AigBig(AjAx)s(BiBx)s — 
A jsB js( Ax Ai) 
— + + 
(123) 


which is to be summed over all possible sets of three terms from the h 


.terms of the original equation (111), allowing repetitions for the reason 
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already noted under m2; that is 2, j,k = 1,2,--, 4. Such a summa- 
tion will hereafter be referred to as summation over the extent of the 
equation,'* and is implied in every symbolic equation. 

In general we may express m, symbolically by the equation 


plmp= ((arb’; | | a’2be,---, a,b’, | a’ (124) 


Expanding by the definition of Art. 5, the leading term is the product 
of p factors of the form a,b’,: a’,b,, because the leading term is made 


without interchanges among the consequents. These factors are the | 


same as @,-a@’,-b,-b’, and symbolically the same as A,gsB,s. Thus if 
we have a set of p subscripts 2, j, k,---, 7, s, t, the t@:ding term in the 
development of p/m, is 


A isBigAjgB ssB ssAisBis (125) 


which is to be summed over the extent of the equation (11). 

The other terms in the expansion of (124) are obtained by making 
all possible interchanges among the consequents, according to the 
definition of Art. 5. Therefore so long as we maintain the polyadic 
notation any term is a product of factors of the form a,b’,:a’,Ds, 
which is the same as a,-a’.,b’,-b,. The vectors a, and a’, might 
correspond to different matrices, likewise b’, and bs. Hence such a 
factor cannot in general be translated into matrix notation if con- 
sidered by itself. For it is of the essence of this transformation that 
every vector a, be associated with its mate a’, and likewise every b, 
with b’,; then a,a’, is symbolically equivalent to A, and b,b’, to B,. 
If, however, by a simple interchange of a pair of consequents in 
the polyadic expression, we obtain a pair of factors (a,b’,:a’.Ds) 
(a,b’,:a’,b,), these may be developed as in (118) and (119) and yield 
the two factors (A,A,;)s(B,B;)s. Such factors occurred already in the 
second term of mz and in the second, third, and fourth terms of m3. 

If three consequents change places among themselves, we obtain a 
product of three factors of the form (a,b’,:a’,b,) (a,b’;: :b:) 
(a,b’,:a’,b,). Itis important to notice that the order of subscripts in 
the vectors a is the same as that in the vectors b but the order of accents 1s 
reversed. It follows that when we develop after the manner of (122) 
the vectors a yield the product of dot products @,- a’ :&1, 
while the vectors b yield the product b’,-b,b’.-b:b’:, b,, where the 
order of subscripts is the same, but the accents are on the first vector in 
each dot product, instead of the second. ‘These factors are equivalent 
in matrix form to (A,A;A,)g and (B,B,B,)s, where the order of sub- 


= 
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scripts in the postfactors B is the same as in the polyadic term, but the 
order of subscripts in the prefactors A 1s the reverse of the order in the 
postfactors B. For in the polyadic expression b’, follows its mate b,, 
and b’; follows b;, and if we rewrite the factor in the form (b,b’,- 
b.b’.-b,b’:)s then b’, also follows its mate b,, whence directly the 
matrix factor (B,B,B,)s._ But in the polyadic expression as it stands, 
each a’, when we follow the same steps, precedes its mate on account 
of the reversal of the order of accents. Thence follows the reversal 
of order of matrices in the corresponding factor. Transformations 
of this type already occurred in the last two terms of m3. 

All the italicized statements in the above discussion are true no 
matter how many consequents have changed places among themselves. 
We may therefore write a general rule for the formation of m,. 

Rule for forming p/m, where m, is the coefficient of 0” in the 
Hamilton-Cayley equation for 6. 

Let there be p subscripts 7, 7, k,---, r, s, t, each of which may have 
any value from 1 to h. Choosing a particular set of values for these 
subscripts, we form a group of p/ terms as follows: the leading term 

is The other terms 
are formed from the leading term by first interchanging the post- 
factors B in all ways, while the prefactors A are at first left fixed in 
position. If a particular postfactor B is left in position, it yields in 
the corresponding term a factor Bg and its prefactor yields Ag pre- 
cisely as in the leading term. If a pair of postfactors as B, and By 
change places, there results in the corresponding term a factor (AA q)s 
(B.B,)s. If three postfactors, as B;, B,, and B,, change places so that 
their new order is B,, B,, By, there results in the corresponding term a 
factor (A;,A,A,)s(B,B.B;)g where the order of prefactors is the reverse 
of the final order of the postfactors. In general if any group of post- 
factors change places among themselves so that their final order is 

there results in the corresponding term a factor (A;A,-+:A,Ax)s 
(B.B,---B,B;)s where the order of prefactors is the reverse of the 
final order of postfactors. 

The group of p! terms thus obtained is to be summed over the extent 
of the linear matrix equation. 


18. CHARACTER OF THE COEFFICIENTS AS ALGEBRAIC POLYNOMIALS. 


It is evident from the form of the leading term that in every term 
of the expansion of m, will be found p prefactors and p postfactors, 
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whether alike or different. Hence m, is homogeneous and of degree 
2p in the scalar elements of the given matrices which define 6. 

Furthermore m, is non-homogeneous and of degree p (in general) in 
the scalar elements of any particular matrix. For if the p subscripts 
which occur in a group of terms are all alike, the p prefactors and the 
p postfactors are all alike. 

Again, each group of terms, corresponding to a particular choice of 
p subscripts, is homogeneous in the scalar elements of each of the 
matrices which occur in the group. The degree in any one matrix is 
determined by the number of times its subscript has been chosen in 
making up p subscripts. The degree in 4, is the same ag that in B,. 

We thus see that each group of terms formed by the above rule is 
made up of terms essentially wnlzke the terms of all the other groups, 
each group, in other words, constitutes a distinct polynomial in the 
scalar elements of certain of the given matrices. Within the group, 
however, various terms may be like one another, and simplifications 
may occur, particularly when special values are assigned to some of 
the matrices; for example, when some of them are allowed to become 
the idemfactor. 


19. THe Equation oF Extent Unity AND ORDER Two. 


As an example, closely allied to that of Art. 10, let us take all given 
matrices to be of the second order, and, to begin with, let h = 1, so 
that the linear matrix equation reduces to the simple form 


AzB =C (126) 


and can be at once solved as x = A~’CB-', Forming the m, by our 
rule we have m, = AgBg which is (a, + dee) (by + be). Next 
2m, = A*sB*s — (A?)s(B?)s. 

The scalar (A?)s may be expressed in terms of simpler quantities by 
using the Hamilton-Cayley equation for the matrix A. Let the 
coefficients in this equation be Ag and A’’s so that 


A? = AgA — A” 
where J is the identical matrix of the second order. Taking the scalar 
of both sides of the equation we have 

= — (128) 


and a similar equation for (B?)s. Substituting in the expression for 
2me we find on simplifying 


ef 
he 
= 
t 
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Ms = — 2A” + A” sB’s (129) 


Since the matrices are of the second order, A’’s and B’’s are the de- 
terminants of their respective matrices; hence on substituting, 


m2 = (ay + (Dirbe. biobe1) 2 12021) (Dirbee 
+ (bi; + boo)? 
+ — (130) 
Again 


6m3 = A®sB%s — 3AgBs(A*)s(B*)s + 2(A5)5(B)g (131) 
by the rule.* But A? = AsA? — A’’gA hence 
(A3)g = — (132) 


by (128), and similarly for B. Substituting values in (131) from (128) 
and (132) to get rid of scalars of powers of matrices we find 


m3 = AsBsA"sB"s (133) 
(ay; + (bu + boo) (11029 A221) by2be1) 
and finally 24m, = A‘sB*s — + 3(A?)2s(B?)?s 


+ 8AsBs(A*)s(B*)s — 6(A*)3(B*)s (134) 
But At = AgA*? — A” gA? hence 
(A*)s = Ats — 4A?g As + (135) 
so that by eliminating the scalars of powers from (134) 
= A”? (140) 
Now in the present example ?C = A°CB?, = A®CB3, etc. Hence 


+ A” + AgBgA*CB? — (141) 
Simplifying by the aid of the Hamilton-Cayley equations for A and B, 


t= 


[AsBsC — BsAC — + ACB| 


checking with the known solution. 


>: 
} 
a, 
= 
- 
‘ 
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20. THe EquaTIon or Extent Two AND ORDER Two. 


The following considerations will suggest how the work of solution 
may be arranged when the equation is of higher extent. We have 
seen that any coefficient m, is a sum of groups of terms. Each group 
is a homogeneous polynomial. Suppose a choice of subscripts in 
which 7 occurs a times, 7 occurs 6 times, etc. Let the development of 
terms, under the rule, corresponding to this choice of subscripts, be 
denoted by G(i%7®---). Each choice of subscripts will occur a number 
of times equal to the coefficient of the corresponding term in the 

expansion of 


(142) 


Thus the entire development may be systematically carried out. 
For example, take 


A xB; + A ;xB; =C (143) 
and form m4. With the above notation we shall have 
= + + + + G(z*) (144) 


By the preceding example we have G(i*) = A” ?sB” 2s and similarly 
for G(j*). By the rule we have 


G(i3}) = A js — 
— 3AjsA js(A?;) isB j3(B?i) + 
+ 6A ;)sBis(B?:B + 24 
— 6 (A*;A;)s(B*:B;)s (145) 


Brevity will be gained in notation, while nothing is lost in explicitness, 
if, in such expressions, we indicate only the matrices A which are 
prefactors in the original equation, remembering that the order of 
subscripts among postfactors, when more than two matrices are 
multiplied, is the reverse of that for prefactors. With this under- 
standing we may also omit the letter A and the subscript S. Thus 
(145) may be abbreviated 


= — 307 (7) (17) — (145 a) 


As a check, if A; and B; are replaced by the idemfactor of the second 
order, should reduce to because mg becomes the 
determinant of A( )B + ( ), and the terms of the third degree in A 


ig 
- 
= 
he 
‘ 
be 
= 
4 
x 
- 
Ry 
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must be the coefficient m3 for the example of the last article. In fact 
we then have, remembering Js= 2, : 


which is the same as (131). 

Developments like (145) are the same in form no matter what the 
order of the matrices involved. They may in general be simplified 
by the use of the Hamilton-Cayley equations for A and B. 

G(ij?) may be obtained from G(2*7) by interchange of subscripts. 
For the middle term of (144) we have 


= — — 44909) + @) (FP) + 209)? + 4007’) 
+ — 2(ajiz) — 4(2°7?) (146) 
which by letting A; and B; be I of the second order should become 


2.A*3B*s — 2(A?)s(B?)s as is easily verified. 
In simplifying, use is to be made of the identity 


(A;A;)s = AigAjs — A*;A; (147) 


where the last term is the star product of A; and A; as already defined. 
Since this scalar is an invariant of the two matrices it may well be 
abbreviated A*;;. Collecting and reducing results we find 


mg = A’? 5B’? + A” is BY jg A* ;B* + A” ;sB™;; 
+ B” 3B" gA%;; — 2A" gB" ig + A" ;gB" jg A* ;;B* 
+ 3B’? (148) 
and by similar processes, 
m3 = AisA” isBigB” ig + AisBjsB” igA*i; + AjsBisA” igB* i; 
— AjgB;gA" isB” is + AjsBjgA” jsB” jg + AjsBigB” jg A* iis 
+ jsB* — AigBisA” jsB" js (149) 
A*;sB” — 2A” isB” is + A” igB?ig + AigAisB* i; — i; 
+ A*jsB” jg — 2A” jsB” jg + A” (150) 
m= AisBis + AjsBjs (151) 
whence (143) is completely solved. 
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1 This symbolic method is due to H. B. Phillips: Some invariants and co- 
variants of ternary collineations,’’ American Journal of Mathematics, 36, 1914. 
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5 For the elementary theory, Bocher’s Algebra may be consulted. 

6 See note 3. 

7 For the laws of p-way determinants see Amer. Journal of Math. 40, 1918, 
by Lepine Hall Rice. 

8 Compare Joly’s Appendix to Hamilton’s Elements of Quaternions. 

9 This result is an extension of Hamilton’s invariant property of his coeffi- 
cients. 

10 Elements of Quaternions, 2nd Ed., Vol. I, Art. 348. 

11 In a series of papers over many years. "They are all, I think, listed in the 
bulletins of the Quaternion Association. See, in particular, Wien. Ber. 112, 
1903, pp. 645, 1091, and 1533. 

12 The term “extent” is due to Sylvester. I have elsewhere given a sketch 
of the present method in its relation to the work of Sylvester. Proc. Nat. 


Acad. of Sci. 8, April, 1922. 
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